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Abstract. The dynamics of excited phonon spectra of the E®e Jahn-Teller (hereafter, JT) model mapped 
onto the generalized Calogero-Moser (gCM) gas of pseudoparticles implies a complex interplay between 
nonlinearity and fluctuations of quasiparticle trajectories. A broad crossover appears in a pseudotime 
(interaction strength) between the initial oscillator region and the nonlinear region of the kink-train lattice 
as a superlattice of the kink-antikink gCM trajectories. The local nonlinear fluctuations, nuclei (droplets) 
of the growing kink phase arise at the crossover, forming a new intermediate droplet "glassy" phase as 
a precursor of the kink phase. The "glassy" phase is related to a broad maximum in the entropy of 
the probability distributions of pseudoparticle accelerations, or level curvatures. The kink-train lattice 
£j ' phase with multiple kink-antikink collisions is stabilised by long-range correlations when approaching 

i— i, a semiclassical limit. A series of bifurcations of nearest-level spacings were recognised as signatures of 

pre-chaotic behaviour at the quantum level in the kink phase. Statistical characteristics can be seen to 
confirm the coexistence within all of the spectra of both regularity and chaoticity to a varying extent 
, (nonuniversality). Regions are observed within which one of the phases is dominant. 

' PACS. 31.30.-i Corrections to electronic structure - 63.22.+m Phonons or vibrational states in low- 

dimensional structures and nanoscale materials - 05. 45. -a Nonlinear dynamics and chaos - 34.10.+X 
• ■ General theories and models of atomic and molecular collisions and interactions 

t-H ■ 

t-h ; 

£^ ■ 1 Introduction as the generalized Calogero-Moser gas (gCM) of Coulomb 
S". \ interacting pseudoparticles 4,5,6,7. This method pro- 
> ' The analysis of irregularities of quantum spectra is es- vides a brid S e to a statistical description of the system and 
K> sentially the search for fingerprints of chaotic behaviour thus appears to be a basis for differentiating between reg- 
^ ■ at the quantum level. It is especially delicate for models ularit y and chaos in some systems [EHTU]. It is worth not- 
ed lacking a reasonable semiclassical limit due to the pres- m S that the statistical description applies only for models 
ence of quantum tunneling or nonadiabatic fluctuations w ^ tn appropriate Gibbs measures, 
[lj. It is equally delicate for systems with a combination 

of regular and chaotic phase space {23 • Previous studies In this paper, the authors seek to numerically and an- 

have shown that nonadiabatic fluctuations act to stabilise alytically examine level dynamics of levels as functions 

the ground state, causing the ground state potential to of the parameter (electron-phonon coupling as a pseudo- 

become highly nonlinear 3 . In this paper we shall inves- time) in pseudospace. Probabilistic considerations enable 

tigate the consequences of the interplay of nonadiabatic determination of a measure of stochasticity within this 

fluctuations due to level correlations and nonlinearity in irregular and highly non-universal model, in the context 

complex excited spectra of the E<§5e JT electron-phonon of earlier analyses of the nearest level spacing probabil- 

model in the situation of increasing interaction strength. ity distributions [TT] . Besides the numerical analysis we 

The JT model, with one control parameter of electron- shall apply an alternative approach of level dynamics. We 

phonon interaction strength, can be mapped onto a fully shall present a related approximate analytical calculation 

integrable classical many-body system. This system is ob- to shed light on spectral structures, the physical mecha- 
served to experience repulsive long-range interactions known nisms of their formation and the crossover between them. 

One advantage of this approach is the ability to classically 

a e-mail: eva.majernikova@savba.sk represent quantum fluctuations and examine nonlinearity 

b e-mail: serge_shp (at) yahoo . com of the originally quantum system. 
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Excited energy levels E n (a), n = 1, . . . , N of a quan- 
tum Hamiltonian in the form H(a) = Hq + aV can be 
considered as dynamic coordinates of iV interacting clas- 
sical pseudo-particles moving in a pseudo-time r = a, 
E n (a) = x n {r). The gCM set of equations equivalent to 
the quantum mechanical problem H(a)\4>) = E n (a)\4>) 
reads [U[5] 



dx /dr-D M — - 2 V LnmL 



nrn-^mn 



dL r , 



dr 



l 



where 



i(r) = ^„(t) - a; m (r)J ■ V mn {r) 



,(1) 



(2) 



p n = (n(r)\V\n(T)) = V nn (r), V nm {r) = (n(r)| V|m(r)>. 

This system is known to be completely integrable [M7J- 
Hence, it possesses as many integrals of motion as it has 
independent variables. The most important two additive 
integrals of motion are the total "energy" (the classical 
Hamiltonian which generates the above equations of mo- 
tion (1)) 



E 



1 N 

2 XX" E ^ 2 fe (r)log|^-x fe | (3) 

n=l l<j<k<N 

and the total "angular momentum" [5 



E 

l<j<k<<N 



(4) 



The first or both integrals of motion were used to develop 
the statistical description of the gCM system by an ap- 
propriate choice of canonical ensemble and corresponding 
Gibbs measures [5][T2]. 

Equation (3) allows us to understand the system as 
a two-dimensional Coulomb plasma of non-local time de- 
pendent pseudocharges V}fc(r). The respective set of equa- 
tions alternative to the equations for L mn is implied by 
(1) and (2) to yield [8] 



dVrr, 



dr 



V m l(T)V ln (T) 



l 



l 



XI 



xi 



-V mn (T)- 



(5) 



The formalism of level dynamics (15) may be applied to 
the excited phonon spectrum of the E<g>e JT Hamiltonian 

H = {b\b 1 + blb 2 + l)I + a{b\+b 1 )a z -a(bl + b 2 )<j x . (6) 

In this equation, &i,2 represent boson (phonon) opera- 
tors of two oscillators with the frequency fl — 1, a x = 



°1 a v = i [ i q J ' a * = U _°! j are Pauli 

matrices, / is the unit matrix, a is the coupling constant 
between the electron and phonon modes. The 2x2 ma- 
trix form accounts for two electron levels. According to 
the above, the model (6) with one non-integrability pa- 
rameter a (H(a) = Hq + aV) and the integrable part Hq 
is suitable for gCM mapping ( 1 ) to be applied. It is worth 
noting that the case containing different frequencies of the 
phonon modes 1 and 2 can easily be transformed into a 
case with two different coupling constants a ^ f3 in the 
two last terms of (6). This was investigated in a series 
of previous papers [3,13,14]. However, the level dynamics 
approach for systems with several nonintegrability param- 
eters is much more complicated [71115). Analysis reveals 
non-trivial gauge structures induced by multidimensional 
parameter space. 

The JT model (6) consists of two degenerate electron 
levels coupled with two vibron modes of different symme- 
try against reflection. It is a typical representation of a 
nonintegrable nonadiabatic system posessing no reason- 
able semiclassical limit. The model has been intensively 
studied over time 3,10, 16 , 17 and is known to be a rich 
source of understanding for physical properties and con- 
sequences, for both heuristic interest and practical ap- 
plications. The rotationally symmetric E(g)e model, be- 
sides of the common SU (2) reflection symmetry, has one 
additional constant of motion - the conserved angular 
momentum J = — 6^62) — \ G y with eigenvalues 

j = 1/2, 3/2, .... Hence, in the following, the representa- 
tion of the definite quantum number j will be used which 
appears to be an additional parameter of the model [llj . 
To illustrate the typical properties of the gCM gas we 
have numerically solved the set of gCM dynamical equa- 
tions (1), (2) and (5) for the system (6). The initial condi- 
tions V mn (0) for the model (6) are specified from the set 
of equations for the excited levels of the model. It consists 
of pairs of even and odd levels E ni n = 2n ri 2n r + 1, 
n r = 0, 1, . . . (main quantum number) determined by the 
well known tridiagonal matrix equations for energy E 



( E 2n r - E)c 2rir + a(/„ r ,„ r c 2 „,. + i + /„ r ,„ r _ 1 C 2llr -l) = 0, (7) 

( E 2n r + 1 ~~ E)c2n r + 1 + a{fn r ,n r C2n r + fn T+1 ,n T C2n T +2) = . 

In this equation, E® represent the energies of unperturbed 
harmonic oscillators; {c„} are components of the wave 
function in the harmonic oscillator representation (Hq is 
diagonal). The perturbations terms in (7) are represented 
by the matrix elements [TT1[T5] 



V2„ + l >2 „(0)= f nrnr = W«r + l + |j-l/2|, 
V2n-l,2n(0) = f nr n r -l = -V^V^V , (8) 

giving the values of pseudocharges (5) at t = 0. 
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2 GCM dynamics of the Jahn-Teller excited 
spectra. Numerical results 

In this section we present a numerical solution for the 
set of gCM equations (12 ) for several values of rotational 

quantum number j. The initial conditions x n (0) = E n °\ 
p n (0) = 0, V2n+i,2n(0), V 2n ~i,2n(0) are given by (8). Here 
and in subsequent sections the index n labels the energy 
levels in increasing order while n r and j represent the 
quantum numbers of the initial unperturbed oscillators 
[TTj . Typically the reduced sets of equations for twenty 
levels were examined. In order to confine the levels we im- 
posed boundary conditions rendering the first and last lev- 
els unmovable, in other words, replacing the first and last 
pseudo-particle by the pseudo-particles of infinite mass. 
The examples of numerical solutions to gCM equations are 
shown in Figures 1 , 2 for three sets of variables (i„ , p n , 
V mn ) involved. The solutions for energy levels not close 
to artificial interval boundaries locally effectively repro- 
duce the numerical solutions for energy levels as functions 
of a obtained via the diagonalisation of the Hamiltonian 
matrix (7) (see [TT] . Fig.l). There the exact level spec- 
trum is illustrated for two different domains introduced 
as "oscillating" and "kink-lattice" domains of the spec- 
trum. These roughly correspond to the different types of 
behaviour shown in Figures 1, 2. As the number of levels 
increases, so does the effectiveness of the solutions at re- 
producing numerical solutions via diagonalisation of the 
Hamiltonian. 

In Figures 1, 2 one can recognise a complex dynam- 
ical behaviour of three coupled groups of excitations: at 
small a there are two groups of pseudo-particles (levels), 
even and odd, of weakly perturbed oscillations with op- 
posite amplitudes (Fig. la) and related rapidly oscillating 
pseudo- momenta p n (Fig. lb). After a small initial period 
the level degenerations are removed due to the creation 
of new inter-level correlations of more distant levels rep- 
resented by the "pseudocharges" V mn (r) for m > n + 1 
(Eqs. (1, 2)). They emerge as the third (central) band in 
Figures 1, 2,c. With increasing a and j the interference 
between the central and side bands increases, destroying 
the oscillatory behaviour of the levels x n and related mo- 
menta p n . This scenario tends to the formation of a single 
stochastic central band of pseudocharges with increasing 
pseudotime (Fig. 2c). Moreover, within the single stochas- 
tic band appear very narrow windows of collision-free bal- 
listic motion, with linear dependence on x n ~ va and 
p n ~ const. This may be attributable to the kink domain 
in the long time limit when the joint effect of the stochastic 
set of correlations (including the long-distance inter-level 
correlations (Fig. 3c)) vanishes, allowing a situation simi- 
lar to free particle dynamics. We shall comment further on 
the structure of the spectra on the basis of the analytical 
results contained in Section |3l 

Correlations between the levels {x n +i, x n ), i — 1,3,5 
(maps) for different model parameters and parts of the 
spectra are illustrated in Figure 3. In these figures the 
notation x n stands for the "reduced" part of the level en- 
ergy with extracted energy of the level Eg and "secular" 



a) j = 63/2 




Fig. 1. Numerical solution of level dynamics equations (1)- 
(2) in the n r > j part of the Jahn-Teller excited spectra for 
j — 63/2 with artificial boundary conditions (see the begin- 
ning of Sect. 2). The levels from n = 93 to n = 102 are shown 
provided that those at n = 92 and n = 103 are kept fixed, (a), 
(b) show respectively "coordinates" x n , and related momenta 
p n = V nn , and (c) shows all the off-diagonal matrix elements 
Vmn of the perturbation V for m, n within the described range. 
The transition between the oscillatory (small a) and the kink 
region with avoidings (p n — 0) is marked by mixing of the Vmn 
branches. The kink lattice region is characterised by strong 
mixing of three V m n branches. Narrow windows of weaker mix- 
ing at large a are apparent. 



part n — a 2 (Cf. Fig. 4 of [TT]). The semi-elliptic shape of 
the maps, related to the highest parts of the spectra, sug- 
gests that the correlations X2 n +i ■ %2n have a strong effect 
with excluded space inside the ellipse (Figs. 3a, 3b). This 
is a direct indication of level avoidings. The effect notice- 
ably changes for different parts of the spectra and becomes 
more pronounced and unpredictable with increasing a and 
j. On the other hand, for extremely high a and j close to 
the classical limit, the effect of correlations randomises 
the map so strongly that the excluded space substantially 
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a) j= 121/2 




b) 




Fig. 2. The same for j = 121/2, n r < j. Note the kink- 
lattice region x n ~ vol at strong mixing of three V nm branches. 
Narrow segments of collisionless ballistic dynamics appear with 
linear time-dependent trajectories x n and constant momenta 
p„ at Vmn = 0. With increasing a the "oscillating" branches 
are absorbed by the central stochastic band. 



shrinks to several regions of a possible fractal structure 
induced by the long-range correlations (Fig. 3c). 



2.1 Probabilistic distributions of pseudo-particle 
characteristics 

A complementary approach to the developed dynamic (de- 
terministic) picture can be provided by describing the 
system in terms of statistical distributions. The common 
ad hoc assumption is to introduce a canonical (or grand 
canonical) ensemble and an appropriate Gibbs measure 
in the phase space of all dynamical variables [51IT2] dM = 
ex.p(—/3H)dx n dp n dL mn . In the domain of developed quan- 
tum chaos, this approach was shown to reproduce the 
main results of random matrix theory (RMT). In partic- 
ular, it is effective in predicting the form of the distri- 
bution of level spacings and level velocities. Similarly, it 
can be used for the investigation of the response of the 
energy spectrum to the change of nonintegrability param- 
eter (that is for the distribution of pseudoparticle veloci- 



ties v = p n = dx n (r) / dr = V nn (r) and their accelerations 
(level curvatures) P(K), K = d x n (r) /dr 2 ). In particular, 
in the domain of applicability of the RMT the distribution 
of the velocities should be Gaussian with the dispersion 
equal to the "temperature" parameter 1/(3 in the Gibbs 
measure j8j. The principal prediction for the distribution 
of the level curvatures is their asymptotic behaviour at 
large K which has been shown universally to be of the form 
P{K) ~ K ~ v with v = 3, 4, 6 for different kinds of ensem- 
bles 12 . The strongest assumption made in this approach 
is determining the form of the Gibbs measure, which im- 
plies that the system is in thermodynamic equilibrium. 
In the domains of undeveloped chaotic behaviour this as- 
sumption appears invalid. In particular, the statistical dis- 
tributions P(v) of the level velocities v = p n = y n „(r) for 
fixed j and r = a for JT system are shown in Figure 4 
They develop from an initial (at small a) oscillation for 
large n (upper part of the spectra, Fig. lc, 4b) to strongly 
broadened oscillatory bands due to the stochastic veloci- 
ties in the lower part of the spectra, Figure 2 c, 4a. Two 
distinct bands of positive and negative velocities and the 
in-gap remnants of the oscillations seen in Figure 4a are 
interpreted in Section 3. 

The non-Gaussian character of the distributions P(v) 
in Figure 4 and P(K), Figure 5 exposes the necessity 
of more correctly introducing the stochastic measure. As 



a) 2=63/2, a=3 



-04 -0.2 



b) j=63/2, a=5 
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c) j-255/2 a=l 
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Fig. 3. Examples of maps of exact solutions x n +i vs x n for i — 
1,3,5 for different j and a and different parts of the spectra: 
levels with numbers 100-300 - blue, 300-500 - red, 500-700 - 
yellow, 700-750 - green. 
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the authors have previously suggested [IT], the develop- 
ment of the distribution of P(K) in "time" a has some 
features common with the diffusion process governed by 
the diffusion equation with additional telegraph term, i.e. 
Td 2 P/dr 2 + dP/dr = Dd 2 P/dK 2 , where D is the diffu- 
sion coefficient. The coefficient of the telegraph term ac- 
counts for memory effects. In particular, it is responsible 
for the appearance of the characteristic "wings" for large 
K (Fig. 5 b) which are not observed through the RMT ap- 
proach. A correct explanation for these features requires 
introducing an equation for which the measure dM(a) de- 
velopes in time and only coincided with the Gibbs measure 
in the limit of RMT ("quantum chaos"). In view of the 
two-peak character of the probability distributions P(K) 
we use their "entropies" S K = - J P(K) log P(K)dK to 
characterize their properties rather than their dispersion 
characteristics [15] . 

The entropies Sk of the probability distributions of the 
accelerations (curvatures) P(K) averaged over successive 
intervals of 100 levels are plotted in Figure 6 as functions 
of pseudotime a. The crossover between the two above 
described phases is evident. These figures show the broad 
maximum irregularity (stochasticity) in the region of the 
maximal mixing of three branches of the pseudocharges 
identified as the region of growing nonlinear fluctuations 
(formation of droplets) as a precursor phase of the kink- 
train lattice (see subsequent section). 



a) a=2j=121/2 n=50-350 

P(v) 

0.07 h 




b) a=2j=121/2 n=350-650 

P(v) 

0.14 



-20 -10 10 20 

Fig. 4. Probability distributions of velocities of pseudopar- 
ticles P(v), v = Vmm(j). Peaked oscillator-dominated region 
in the upper part of the spectra (b); Stochasticity-dominated 
broad bands in the lower part of the spectra (a). 



3 GCM dynamics of the Jahn- Teller excited 
spectra. Approximate analytical approach 

The simplest analytical approach based on the approx- 
imation of constant (initial) values of matrix elements 
V 2n +i, 2 n(i~) and Vm~\,2n{j) can be applied for a small 



a. It must be borne in mind that this approach neglects 
the dynamics of V rnn (r) , the correlations of more distant 
levels developed in the course of pseudotime (Figs.3,b,c). 
As a consequence our two-level approach includes only the 
dominant nearest neighbour interactions. To a limited ex- 
tent this approximation allows for understanding the main 
features of the model and is satisfactorily consistent with 
our numerical results. This formulation recovers the local 
crossover between the anharmonic oscillator region and 
the complex nonlinear one. 

Within the approximation of constant matrix elements 
in (2) it may be taken that V 2 n+i,2n(r) ~ V 2n +i,2n(0) = 

fnn{n r ,j), V 2n -l,2n{T) « V 2 n-l,2n(0) = /n,n-lOv, j) from 

(8) (in what follows we denote j = \j — 1/2] for brevity). 
Then, 

^2«+1.2«(t) = V2n+l,2n(0)(a;2n+l(T) - X 2 „(t)) , 
Lln-l,2n{T) = V2„-l,2n(0)(x2„(r) - X 2 „-i(t)) . (9) 

Since the number of phonons is not conserved the num- 
ber of pseudoparticles (excited levels) in the system with 
Hamiltonian (6) is not constant. Therefore a two-level ap- 
proximation of the gCM equations can be used, but only 
for small a. Denote the corresponding levels as x\, x 2 . 
The origin may be set so that two levels x\ 1 x 2 are sym- 
metric and separated by x 2 — x\ = 2x; the effect of the 
other levels of the system can be accounted for by set- 
ting two bounding (fixed) levels, xq, X3 separated by a 

a) ff=2j=121/2 n=50-350 

P(K) 




-10000 -5000 5000 



b) <r=2j=121/2 n=650-950 

P(K) 




Fig. 5. Probability distributions P(K) for a = 2, j = 121/2, 
where K = d 2 x n /dr 2 is the level curvature (pseudoparticle ac- 
celeration) . Note the two-peaked oscillator-dominated region in 
the upper part of the spectra (b). The stochasticity-dominated 
region in the lower part of the spectra exhibits (i) a prominent 
peak at K = which corresponds to the kinks (p — const, see 
Fig. 2d); (ii) characteristic for RMT prediction of asymptotic 
power- like dependence at large K. The side "wings" of the dis- 
tributions are in contrast to the predictions of the RMT and 
indicate a presence of memory terms in the equation for the 
measure dM(a). 
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S K j=121/2 S K j=63/2 




Fig. 6. Entropies Sk of the probability distributions of the accelerations (curvatures) P(K) averaged over successive intervals 
of 100 levels. The maximum irregularity (stochasticity) appears in the region with maximum mixing of three branches of the 
pseudocharges identified as the region of droplet formation and as a precursor phase of the kink-train lattice. Triangles - levels 
50-150; stars: levels 150-250; squares: 250-350; open triangles: 350-450. The entropies are given with an arbitrary constant of 
the form AK log AK accounting for the finite sampling bins AK on the histograms such as shown in Figure 5 , which were the 
same (AK — 10) for all graphs. 



distance £ 3 — xq 

fnn — /oi fn+l.n 

yields 



3. To simplify notations we denote 
fx. The set of gCM equations then 



f 1 dpi 
2~dr 

1 dp 2 

2 dr 



-/o 2 



fi 



x 2 - X\ 

fl 

X 2 — X\ 



Xi - Xq 

fl 
X3 - X 2 



(10) 



The equation for the separation x of the two nearest 
trajectories reads 



dV[x) 
dx 



2x 



Vl 



3/2 - x 



(11) 



where V(x) = —f§logx — 2ff log(1.5 — x) is the one- 
dimcnsional time-dependent Coulomb potential with x eq = 
3/q/2 (/o + Vi) as a point of the potential minimum 
where V (x eq ) = 0. 

The harmonic oscillations of the frequency V" = 2(2 f 2 + 
/o) 3 /(3/oA) 2 become unstable at V = 0, or when f§ = 
2f 2 or n r ~ j. This condition determines the border be- 
tween the harmonic oscillation region (n r > j) and the 
nonlinear region (n r < j) in this simplified model. The 
model generates the map illustrated in Figure 7, Since 
the subsequent clusters in it are independent, the whole 
map is completely determined by the phase shift of the 
oscillations from cluster to cluster (changing the param- 
eter n r numbering the clusters). Resemblance to the cor- 
responding exact results, Figure 3a can be traced only 
for higher levels otherwise the simple two-level symmetric 
model completely ignores the level avoidings. 

In order to proceed beyond the two-level harmonic 
approximation we shall account for the collective effect 



of neighbour levels with level-to-level variations of the 
"pseudocharges" f m n- Let us define small fluctuations 5 
as X2n+i — X2n = 1 + ^2n- The set of equations derived 
from (7) then reads 



1 d 2 S 2n 


Vnn 


f 2 

J n+l,n 


f 2 

J n.n— 1 


2 dr 2 


l + S 2 n 


1 + S 2n +i 


1 + &2n-l 


1 d 2 5 2n+ i 


Vn+l,n 


f 2 


f 2 


2 dr 2 


1 + $2n+l 


1 + S 2n +2 


i + s 2n 



(12) 



For sufficiently large n, f n +i, n +i ~ fnn, fn+i,n ~ fn,n-i, 
hence equations ( 12 ) imply the approximate relations <52n+i 
~ ^2n-i and ^2ra+2 ~ $2n- Then the pairs of pseudoparti- 
cles (levels) represent an analogy of a dimerized system or 
a lattice with two atoms in a cell. It is worth noting that 
from (12) it follows that 52m$2n+i is a P a h of reflection- 
symmetric variables against the transformation r — > it. 



j=63/2 or=3 







\ ** y 1"..^" * a* 

' • r ^ 




-03 -0.1 -0.1 




-Q.J 




-0.1 




-0.3 





Fig. 7. Mapping of approximate solutions x n +i vs x n (as in 
Fig. 3a) generated by the simple two-level model (10). The cor- 
respondence with the exact result in Figure 3 is satisfactory 
except in the near-zero part since the level avoidings are ne- 
glected in this model. 
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3.1 Structure of the nonlinear region and quantum 
phase transition-like nature of the crossover 

For increasing r and j the picture of "pseudocharges" 
V mn (T) (Eq. (3) and Fig.l) dramatically changes by the 
creation of new long-range components which form a cen- 
tral branch of highly irregular oscillations. A strong mix- 
ing between the short-range (two side-bands in Fig.l) and 
long-range V rnn (central band) proceeds towards a limit of 
large j, when the side-bands disappear (Fig. 2:). At large 
a, the level dynamics is mostly affected by the central 
irregular branches of the "pseudocharges" V mn . A sim- 
ple approximate dynamic description of the crossover be- 
tween oscillatory and nonlinear regions can be obtained 
from the dynamic equations (12) for the fluctuations, 5 2 n 
and #2n+i- We develop the terms on r.h.s. of equations 
(12) into a series up to the lowest order of nonlincarity S 3 . 
In the following approximate calculations we shall move to 
a continuum approximation in the coordinate n labeling 
the level system. If we extract the term 



1 



<^2n+2 + £>2r. 



25 



a 2 * 



2n+l 



2n+l 



dn 2 



J 2n+1 



(13) 



and define <52n+i = 1 + y> S 2n = x, the equation for #2n+i 
of ( 12) can be rewritten in a suitable form for further 
analysis, 



_ " _ f 2 " — f 2 - f 2 _i_ 

2?/ Jn+l.n+lV — Jn+l,n+l Jnn + 

2 (/„ 2 +i,„+i - Vl + i tn ) V - 2fn + l,nV 2 (y + 2) (14) 

— (/n+1, n+1 — fnn) x ~ (/n+1, n+1 + fnn) ( x ~ X )i 

where the index n of y n has been omitted for simplicity. 
Using a periodicity condition mod(y + 2) = y the equation 
(14) is rewritten as follows: 



l.._ f2 »_ dV(y) 

2^ Jn+ln+lV — Qy 

where the nonlinear potential 



1 



(15) 



j < n r , fl +hn+l - Vl+\,n < 0> V(y) = -f%+i <n y - 

(/n+1 n+1 ~ 2fn+i.n)y 2 becomes a single-well potential, 
and we identify the phase of anharmonic oscillators. At 
j > n r one has fn+\n+i — 2f 2 +ln > 0. Here a double 
potential well of V(y) opens, and correspondingly a new 
phase appears. Equation (15) with the driving force omit- 
ted acts as an analogy to a quantum second order phase 
transition for each level (trajectory) n. However, the small 
driving force in (16) changes this scenario before it van- 
ishes at large j and a = r. Examining equation (14) with 
accounting for a constant (y-independent) driving force 

F = [-(fn+l,n+l - fin) + (/n+l,„+l + fLV] (!"*) > 

(17) 

results in 



2^ fn+l.n+lV ^(fn+l,n+l 2/^ + i iTl )2/+2/ 2 +1 „J/ 3 — F, 

here, F in equation (17) consists of two competing parts: 
the drift term due to the broken periodicity of the level 
spectra at finite r and the "pressure" force from the neigh- 



bouring levels. If F > and f 2 +1 n+1 



2/S+i lt , > h 



drives a first order quantum phase transition-like crossover 
generating nonlinear fluctuations. The equation of type 
(18) in the dimension d > 1 has been investigated [19 
as a prototype for a first order phase transition involving 
the coexistence of two phases, one phase represented by 
growing droplets (bubbles) spanned by the force F in the 
sea of the other competing phase. In our case the dimen- 
sion d = 1; Nevertheless, there exists an exact solution 
to the normalized nonlinear equation of the type (18). It 
was found in the closed form of a nonlinear oscillation [20 j 
which in our case reads 



v ' w ' 2r ~ 7j ' ' = a 



ni + cos ( w ( 2r 



ri9 + cos [ w 2r = 



where 



(19) 



v(y) = 2fl +ln y A - {fl+i, n+ i - Vl + i, n ) y 2 - (16) 

(/n+l,ra+l — fnn) (1 — X )V + (/n+1, n+1 + fnn) ^ (1 — x ) V 

indicates a crossover analogous to the quantum first order 
phase transition due to a small linear driving force ~ y, 
(\x\ <C \y\ by definition). With the simplification valid for 
large n, f 2 n w f? l+ln+1 , fl +ln « / 2 „~i one can obtain a 
condition defining the dimerisation region &i n w — ^2n-i 
where equation (12) for &i n yields an equation analogous 
to (15). In the case when the driving forces on the r.h.s. of 
( 16) (~ y) are either mutually compensated or both tend 
to zero (large j and n r ), that is, (fn+i, n +i - fnn) { l ~ x ) ~ 

(/n+i,n+i + fnn) x2 {^~ x ) < /™, equations (15) and (16) 
imply a crossover analogous to the quantum second order 
phase transition at / 2 +ljIl+1 - 2f 2 +1 n = 2(j - n r ) = 0. 
This coincides with the criterion for the crossover to non- 
linear region from simple analysis of equation (11). For 



"l 



2 ^ n 2 2 _ /n+1, n+1 ^fn+l,n n 



71 2 



-, a 



f 2 2 + n 2 > ° ' 



= (fn +l,n+l 



II 



/n+1, 



n+1 



2/n+l, 



2 

> 0. 



The "droplet" fluctuation ( 19) is a traveling periodic 
non-sinusoidal oscillation moving with velocity v ~ 
The amplitude a of the fluctuation y is spanned by the 
driving field F. This force grows when x (x < x max = 
1 / y/2(j + n r )) decreases until the energy of the excitation 
(19) reaches the energy of the kink, where v w 2-y/J is 
the minimum velocity of the kink. If j increases, j S> n r , 
then x — >• and a — > (a — > l/(j — n r )), so that the 
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intermediate "droplet" glassy phase transforms into the 
more stable kink lattice phase (see below). 

Simultaneously, with increasing j and a the dimeriza- 
tion condition mod(y + 2) = y fails (Fig. 2), f 2 +1 n+1 ~ 
/„„ and the driving force F vanishes (x max — > 0). The 
equation (14) up to second order terms tends to 



2^ fn+l,n+lU — 2(/„ +1 n+1 2f n +l,n)y ^fn+l,nV ' 

(20) 

If we move to imaginary space C = i(vr — n), then the 
solution of (20) yields a traveling pulse 

^ n+1 — Vn+l.ii) 



4/n^+l,, 



x cosh 



/n+1 



+ l,n+l 



2 /„+!,. 



2/n+i 



n+1 



1/2 



(C-Co) .(21) 



From the definition of $2n+i we have X2n+2 — ^2n+i ~ 
dx2n+i/dn = 1 + <52n+i = 2 + y and equation (21) implies 
a kink-shaped profile, 



x 2n +i = const + 2n 
3 




(/n+1, n+1 2/ n+l! „)(2/ rl+l n+1 V ) 
1 1/2 



1/2 



/n+1, n+1 Vn+l,n 
2 /n+1, n+1 — l ' 2 



(C-Co) , (22) 



where v 2 < 2f 2 +ln+1 
clling kink profile and 



is the velocity of the trav- 



2/ 



n+1, n+1 



/n+1, n+1 



+ l,n 



1/2 



is the kink width. The time of collision can be estimated 
fr om vt c = 2L> (2f 2 + hn+1 / (f 2 +hn+1 - 2/„ 2 +1 ,„)) 1/2 ~ 

\J (2(j + n r ) j (j — n r )), j > n r . Hence, the collision time 
close to the semiclassical limit yields r c — * 1 / \/2j. 

Equation (15) with omitted linear terms in the poten- 
tial (16) represents a nonintegrable member of a family 
of Klein-Gordon equations, that is, <£ 4 equation known 
to exhibit kink, antikink and a double kink solitary solu- 
tions [23]. The solution (22) represents the tunneling be- 
tween two nearest neighbor trajectories. Indeed, for large 
3 (/n,n > 2 /n+i,n or 3 > n r) it gets a form of a prop- 
agating nonlinear solitary excitation, a kink (antikink), 
in the imaginary space y(£ — Qq), £ = i(vr — n), where 
Co = i(vTo — uq) restores the translation symmetry. 

The transition to the imaginary space effectively re- 
verts the sign of the potential, that is, a pseudoparti- 
cle turns to the tunneling domain and changes its parity 
(cf. the mirror symmetry of equations (12)). For example, 
<Wfi(ui) -> 8 2n (vi), S 2n (v 2 ) 8 2n +i{v2)- As a result, 
the trajectories (pseudoparticles) interchange their veloc- 
ities vi,v 2 : S 2n {v 2 ) -> S 2n (vi) and S 2n +i{vi) -> S 2n+1 (v 2 ) 



when transferring from the real to the imaginary space 
and vice versa. This scenario can be understood as a se- 
ries of kink- antikink collisions related to two subsequent 
levels. 

Thus, the nonlinear phase represents the gCM trajec- 
tories as a kink-antikink chain structure apparent in Fig- 
ure 2. Such a solution qualitatively agrees with the re- 
sults of authors dealing with similar systems within the 
gCM level dynamics approach 6,8,21,22] predicting the 
"soliton"-like structures. 

The real configuration of the level system under con- 
sideration is described by the discrete label n. In recent 
years more realistic discrete lattice versions (13) of Klein- 
Gordon equations have been in the focus of interest; they 
show a rich variety of behaviour. Besides the resonance 
energy exchange at kink-antikink collisions, a recently dis- 
covered aspect of this behaviour is the chaotic scattering 
upon certain conditions and fractal behaviour of veloci- 
ties after collision strongly dependent on initial conditions 
24, 25, 26, 27, 28 . The kink-antikink collisions in our sys- 
tem might follow this scenario of initiating the chaotic 
behaviour and fractal structure at large j, Figure 3c. 

In the case of d 2 S/dr 2 = we have p n = const so that 
the value of the curvature K n = d 2 x n /dr 2 — which 
corresponds to the peak in Figure 5 a represents the kinks. 
One has x 2n +\ = x 2n + vAt (here At is the time distance 
between two kinks) and from (12) one finds 



>2n+l 



2/rm + / 



f 2 

J n,n—l 

I 2 

Jn+l,n 

2 i j'J 

n+l,n 



<^2n-l 



fnn—1 2 t /*, 



2 

nn 



/n+1, 



f 2 

J n+ln 



<5 2 n(l 



(23) 



and 



o f 2 _ f2 _ f2 

Jn+l,n Jn+l,n+l J nn 



+ 



t-J. ~ 2f 2 

./n+1, n+1 x MX \ i Jnn 



(24) 



2/n+ln 



'^2n+2(l — &2n+2) + 



2/n+li 



-#2n(l - S 2n ) . 



Equations (23) and (24) are reminiscent of the logistic 
equation x n+ i = Ax n (l — x n ) which implies a known 
scenario for the transition to the chaotic region for A = 
A crit = 3.56994 ... In our case A ~ (f 2 +hn+1 + f 2 n ) /Vn+i.n 
~ (n r + \j — 1/2|) /n r . The ranges of values of j and n r 
considered above give the values of A > 2 close to 2 from 
above, that is, from the subcritical region. The hallmark 
of chaos at the quantum level can be identified as a pre- 
chaotic region of the series of bifurcations of the interme- 
diate glassy phase. 

The pre-chaotic behaviour applies to the medium re- 
gion of values of j and n r between the domains of weak 
coupling with j < n r (dimerized pairs of damped oscilla- 
tors) and strong coupling with j 3> n r . The description 
via initial quantum numbers at a = 0, n r and j related 
to the radial and rotational degrees of freedom of the un- 
perturbed oscillators, is kept by the mapping onto the 
classical gCM gas. For large j, j ^> n r , the short-range 



E. Majermkova, S. Shpyrko: Level dynamics approach to the Jahn- Teller model 



9 



correlations are strongly weakened while the long-range 
correlations dominate, as well as fluctuating around zero 
(Fig. 2c). 

The kink-antikink chain level structure describes the 
special kind of self-avoiding motion of the gCM trajec- 
tories in the integrability parameter a as a pseudotime. 
Alternatively, one can also consider the levels as moving 
along almost straight lines broken by small intervals where 
they undergo "collisions" with neighbouring levels. They 
can be clearly visualized as a series of propagating pulses 
through the series of gCM pseudoparticles, along straight 
lines of positive and negative velocities. The apparent reg- 
ularity of this structure allows us to introduce the idea of a 
superlattice of excitations in the space of gCM pseudopar- 
ticles. It is this new structure which we call the kink-train 
lattice, which may represent a new quantum chaotic phase 
in certain segments of the JT spectra, created by the kink- 
antikink collision mechanism described above. 

A multiplicity of kinks and antikinks (solitary waves) 
travelling through gCM particles may be readily observed 
in Figure lb of the another paper on the subject by these 
authors [TT] . It is recognisable in the numerical gCM sim- 
ulations in Figure 2 a, and in the statistical distribution of 
velocities of gCM pseudoparticles, Figure 4. The lack of 
symmetry of the positive and negative velocities against 
zero illustrates the driving of the lattice by a field su- 
perimposed by the broken periodicity, in turn implied by 
/n+i,n+i > fnn over the whole lattice. 

The trajectories forming the kink-antikink train lattice 
move as almost free particles between the collisions, with 
velocities of opposite sign distributed within two bands 
(Fig. 4a). The remnants of oscillations within the gap 
demonstrate the coexistence of both oscillating and kink 
phases. The long-range order developed at large a sta- 
bilises the almost periodic kink-train lattice (Fig. 6). 

It is worth noting that for j > n r the effective val- 
ues of the logistic map constant A can reach the values 
above the critical A c for some n r . This indicates chaotic 
behaviour of the level spacings. However, the cumulative 
spectra (including all j) impose intersections of the trajec- 
tories because of the multiple degeneracy in j which will 
dominate and imply, for example, the Poisson-like prob- 
ability behaviour of the nearest neighbour level spacing 
distribution P(s) at s — ► [13] • Hence we can conclude 
that a coexistence of regular and irregular nonuniversal 
behaviour is characteristic for the excited spectra of the 
present model and in the whole extent of the interaction 
parameter. 

3.2 Semiclassical approximation 

One possible way of identifying the quantum chaoticity 
is to provide an analogy to the behaviour of the system 
at a classical or semiclassical level. The peculiarity of the 
class of electron-phonon systems is the marked ambiguity 
of passing to a semiclassical approximation from an initial 
set of Heisenberg equations implied by Hamiltonian (6). 
The common algorithm, decoupling products of quantum 
observables and replacing them by the classical ones, can 



be performed in different ways. The problem of moving to 
a semiclassical approximation was touched in a previous 
paper [11]. Nevertheless, the investigation of trajectories 
of semiclassical observables reveals a picture of classically 
chaotic domains. These at least present a qualitative cor- 
relation to the picture of "quantum chaotic" behaviour 
outlined above. Respective equations of motion are im- 
plied by the semiclassical version of the Hamiltonian (6) 
in the form 

H = ^(P 2 + Pf + Q\ + Ql) + aQ lZ - aQ 2 x . (25) 

Here, we introduced the coordinates of phonons Qi, 
Pi and of the electron x = (a x ), y = (<r y ), z = (<x 2 ). 
The electron coordinates satisfy the condition x 2 + y 2 + 
z 2 — 1 for the Bloch sphere. Additionally, besides the 
conservation of the energy (25), the angular momentum 
J = Q1P2 — QiP\ — y is a constant of motion. As a result, 
the system has four degrees of freedom. 

For illustrative purposes, in Figure 8 we present the 
Fourier spectrum of the semiclassical trajectory Q\{t) (a 
picture for Q2 looks similar) for different values of the 
energy £ of a system. The middle part of the energy in- 
terval (corresponding approximately to the intermediate 
part of the quantum spectrum with the strongest mix- 
ing of three effective potential wells) shows the marked 
domain of classical chaos. Meanwhile, the motion of the 
classical trajectories at low and high energies is noticeably 
regular. 



4 Conclusion 

The dynamics of the excited levels of the E(gie Jahn- Teller 
model, mapped onto the generalized Calogero-Moser gas 
of interacting pseudoparticles and developing in pseudo- 
time, shows a complex interplay between the nonlinear- 
ity and fluctuations of the dynamical degrees of freedom. 
Numerical results of Section 2 on lattice dynamics and 
the probabilistic characteristics were interpreted, at least 
qualitatively, on the basis of the approximate analytical 
results of Section 3. 

In the intermediate range, j ~ n r , the gCM dynamics 
implies a maximum degree of irregularity caused by the 
interference of three branches of the pseudocharges with 
growing range of interaction. We have found that the max- 
imum mixing, and hence the broad maximum of stochas- 
ticity (irregularity) illustrated by the entropy of curva- 
tures of the gCM trajectories (accelerations), is related 
to the new intermediate glassy phase of nucleation. This 
is thought to represent the growing kink-antikink chain 
structure. This scenario can be thought of for each level as 
an analogous to a series of local first order quantum phase 
transitions. In the intermediate droplet glassy phase, the 
phase of damped (anharmonic) oscillations dominating at 
j < n r and the lattice of interacting kink-antikink chains 
of the gCM trajectories dominating for j > n r coexist. 
As an alternative picture of this chain structure, we have 
recognised the kink-train lattice phase of the levels. These 
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Fig. 8. Frequency spectrum of the Qi trajectory in semiclassical approximation for different values of total energy 



move along straight lines broken by small intervals, where 
they undergo "collisions" (avoidings) with neighbour lev- 
els. They form a superlattice of excitations in the space 
of gCM pseudoparticles which we believe to represent a 
new quantum chaotic phase. The trajectories forming the 
kink train lattice of the kinks and antikinks move almost 
as free particles between collisions, with velocities of op- 
posite sign distributed within two bands (Fig. 4 a). The 
long-range order developed at large a stabilises the al- 
most periodic kink-train lattice (Fig. 6). It is worth re- 
membering that throughout this work we used the rota- 
tional quantum number j as a parameter. Real spectral 
characteristics are cumulative ones including the contri- 
butions of all j's, where the multiple degeneracy in j (level 
crossings) brings in components of regularity. Hence, real 
statistical properties present the mixing (coexistence) of 
contributions with prevailing higher degree of regularity or 
of those with more irregular nature. These are, however, 
highly nonuniversal. 
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